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Abstract Nomenclature 


The current operational practice in scheduling air 
1 raffie arriving at an airport is to adjust flight sched- 
ules by delay , i.e. a postponement of an aircraft’s 
arrival at a scheduled location, to manage safely the 
FAA-mandated separation constraints between air- 
craft. To meet the observed and forecast growth 
in traffic demand, however, t fie practice of time ad- 
vance (speeding up an aircraft toward a scheduled 
location) is envisioned for future operations as a 
practice additional to delay. Time advance has two 
potential advantages. The first is l lie capability to 
minimize, or at least reduce, the excess separation 
(the distances between pairs of aircraft immediately 
in-trail) and thereby to increase the throughput of 
the arriving traffic. The second is to reduce the 
total traffic delay when the traffic sample is below 
saturation density. A cost associated with time ad- 
vance is the fuel expenditure required by an aircraft 
to speed up. We present an optimal control model 
of air traffic arriving in a terminal area and solve it 
using the Pontryagin Maximum Principle. The ad- 
missible controls allow time advance, as well as de- 
lay, some of the way. The cost function reflects the 
trade-off between minimizing two competing objec- 
tives: excess separation (negatively correlated with 
throughput) and fuel burn. A number of instances 
are solved using three different methods, to demon- 
strate consistency of solutions. 
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aircraft 

the number of flights 

the index of the a/c on the i-th flight 

the aircraft type of the a/c on the i-th flight 

the arc length coordinate of the z-th a/c 

the vector (si, S 2 , . . . , s/) 

the minimal cruising speed of the z-th a/c 

the maximal cruising speed of the z-th a/c 

the actual or assigned speed of the z-th a/c 

the vector (v\, V 2 , ...,?;/) of control speeds 

fuel burn for an a/c of type r(z) at speed v 

the “economical cruising speed” of a/c z, at 

which the fuel burn is minimal 

exit time for problem PI (definition 2) 

constant parameters in B r ^(vi) (see (11)) 

the Euclidean inner product [1, section 5.2-6] 

the norm [1] arising ffom( , ) 

“is defined as” 

required min. separation for aircraft z di- 
rectly followed by aircraft i 
excess separation (13) 

rela 1 i ve importance coelT. of B r ^)(vi ) to E( s) 
the optimal control problem for the Delay- 
Time Advance flight phase (definition 2) 
distance the z-th a/c is to fly from the end of 
the Delay-Time Advance phase to landing 


Introduction 

The current operational practice in scheduling air 
traffic arriving at an airport is to adjust flight .sched- 
ules by delay , i.e. a postponement of an aircraft’s 
arrival at a scheduled location, to manage safely the 
separation constraints between aircraft. This prac- 
tice is most effective when the system is near satu- 
ration or capacity, leaving delay as the only control 
option. Emphasis is placed on delay distribution 
between the Center and the TRACON (see refer- 



ence [2, section 2.2.5] and references therein). In 
this practice, delay is preferred to a speed-up known 
as time advance. The practice of time advance, al- 
though not in use today, is envisioned for future op- 
erations in light of the observed and forecast growth 
in traffic demand. The ability to practice time ad- 
vance can reduce overall system delay in two types 
of situations: i) in the face of impending saturation, 
and ii) when the state of traffic alternates between 
saturation and under-saturation. A cost incurred by 
an aircraft exercising time advance is increased fuel 
expenditure. This calls for a study of the trade-off 
between time advance and fuel burn. 

The focus of this paper is such a study, carried 
out using optimal control theory [3] , and leading to 
a computationally efficient procedure for produc- 
ing speed advisories. Air traffic arriving into a ter- 
minal airspace, with multiple runways allowed, is 
modeled using optimal control, with the state space 
constructed as in [4]; the portion of this modeling 
framework used in this paper is formulated in the 
section A mathematical setting for optimal control 
models with explicit delays). 

A number of past research efforts have focused 
on the trade-off between delay and fuel burn [5,6]. 
An approach widely used in general Air Traffic- 
Management (ATM) research is to cast the prob- 
lem as a mixed-integer (non) linear program [7-10]. 
Since the sets of flights, waypoints, meter fixes, 
and route segments are finite, it is natural to think 
of the operational problems intuitively as discrete. 
The mixed-integer programming framework corre- 
sponds directly to this intuition, hence is a conve- 
nient model for capturing the problem realistically. 
In the context of modeling with a view toward fu- 
ture operational use, however, this framework faces 
the challenges of i) a lack of qualitative insight into 
the behavior of optimal solutions, and ii) absence 
of proofs of convergence or of bounds on the com- 
putational cost, or of both. 

One advantage of optimal control over other ar- 
eas of optimization is the availability of a general 
theory, notably the Bellman Optimality Principle 
and the Pontryagin Maximum Principle [3]. This 
t heory offers insight into the qualitative behavior of 
optimal solutions, often before a solution is com- 
puted, and can yield computational solution proce- 
dures that come with mathematical proofs of cor- 


rectness (i.e., that a computed result is in fact a 
solution to the problem) and of low computational 
complexity bounds (i.e., that the computation will 
complete in a low-degree polynomial time) [11,12]. 
Such proofs are a highly desired feature for an algo- 
rithm that drives an automated tool for supporting 
safety-critical real-time ATM operations. 

The solutions found in this paper for problem 
PI (see definition 2) make the following three main 
contributions: 

• a demonstration that the model is scalable; 

• an efficient way (eligible for real-time use) to 

compute optimal speed advisories accurately; 

• insight into the relation between fuel burn, 
incurred by time advance, and throughput, 
characterized by excess separation in the 
Delay-Time Advance phase of the flight . 

Background 

In general, the arrival of air traffic can be divided 

into the following phases. 

• Delay-Time Advance phase: All aircraft are 
still able to exercise delay or time advance. 

• Mixed phase: At least one aircraft can no 
longer exercise delay or time advance, and at 
least one other aircraft still can. 

• Separation Minimization phase: Aircraft can 
no longer exercise delay or time advance. 

The focus of this paper is on the Delay-Time Ad- 
vance phase. This phase is modeled, below, as an 
optimal control problem which includes a running 
cost term reflecting fuel burn. Two assumptions 
central to the model are as follows: 

Assumption 1 In the Delay-Time Advance phase, 
an optimal trajectory does not touch the boundary 
of the set of conflicting states (i.e., states in which 
at least one pair of aircraft violates the separation 
constraint). 

Assumption 2 The aircraft in question are flying 
in zero-wind conditions. 



One may question whether the mathematical model 
developed here with the use of the simplifying As- 
sumption 2 is readily generalized to include reason- 
able uncertainties, such as wind. The role of As- 
sumption 2 in this study is as follows. Since an 
aircraft’s navigation system will generally keep the 
aircraft on the prescribed path, the main effect of a 
wind field is from its component tangential to t lie 
path. The effect of uncertainties, such as wind or 
error in control execution, on the model is an in- 
troduction of error terms into the endpoints of the 
aircraft’s speed range: 

(min. speed) -1- (error due to uncertainty) 

< (the aircraft’s ground speed) 

< (max. speed) + (error due to uncertainty) 

Trajectories allowed to touch conflict boundaries 
before reaching minimal separation are beyond the 
scope of this paper. Optimal control models cor- 
responding to the Mixed phase and the Separation 
Minimization phase are given, respectively, in ap- 
pendices I and II. 

All the operational decisions concerning the use 
of delay or time advance are made in the flight phase 
corresponding to optimal control problem PI, for- 
mulated below, with a weighted objective function 
of the form 

(coeff. of relative utility) (collective fuel burn) 

+ (excess separation) 

A solution to the problem with a total of I air- 
craft is a time-dependent vector-valued function 

v(0 = (»i(0 «/(«)). (!) 

where Vi(t) is the speed dictated to aircraft i at time 
t. 

Remark 1 Operationally, it is desirable that the 
control strategy v(t) be piecewise constant: in this 
case, every time interval of constancy can be re- 
garded as a speed advisory, and v(t) is then a se- 
quence of speed advisories for the entire set of air- 
craft. 

A general solution is found below as a root of an ex- 
actly known polynomial. For the cost functions con- 
sidered, the mathematical model indicates piecewise 
constant control strategies v(£). Two by-products 
of the solution are as follows: 


• a functional interdependence between the 
speed advisories and the above coefficient of 
relative utility, and 

• a way to compute-once the solutions for other 
flight phases are found the desired change (a 
delay if positive, a time advance if negative) to 
the scheduled arrival time (STA); see the sec- 
tion Obtaining delays and time advances from 
a solution. 

Solutions to a number of problem PI instances are 
presented in a later section. 

A mathematical setting for optimal 
control models with explicit delays 

The state space 

Consider a collection of P paths in R*, indexed by 
the set P = {1, . . . , P}. Path p € V is assumed to 
admit a continuous and piecewise continuously dif- 
ferentiable one-to-one parameterization by its arc 
length s p . Thus, if x p is the mentioned spatial pa- 
rameterization of path p , then 

x p (sp) € R\ — oo < s p < 0 for all p € V (2) 

Definition 1 The point at s p = 0 in parameteriza- 
tion (2) will be called the distal point of path p. 

A pair of paths can either be disjoint, or intersect 
at a point (Figure 1A), or overlap (Figure IB). The 
distal points in Figure 1 are shown as black dots; 
the parameters a decrease, in the direction of air- 
craft motion, toward the distal point. The portions 
of the paths near the crossing (overlap) play a role 
in the geometry of the conflicting states, as will be 
shown in Figure 2. 

Suppose each aircraft in a set A = {1,2 ,...,/} 
is assigned one of the paths in V. To specify that 
a path p £ V has aircraft i assigned to it, we will 
write p(i) instead of p. The following definition of 
a state space follows the approach first introduced 
in [4]. Since the parameterization mapping x p is 
one-to-one and onto, the spatial position of aircraft i 
is completely determined by specifying its distance, 
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Figure 1: Path pairs. 


measured along the path, from the distal point (defi- 
nition 2). Let —Si be aircraft i’s distance, measured 
along the path, from the distal point. This defini- 
tion implies that if aircraft i is on path p(i ), then 
s i — s p(i)- Thus, if each aircraft has already en- 
tered its path, then by specifying all the Si s (each 
Si increases with time), one specifies each aircraft's 
spatial position. Accordingly, in the control model 
constructed below, the vector 

s = i s i)i£A 

is a state vector, hence is regarded as a function of 
time: 

S = s (i) = ( Si(t)) ieA 

It follows that the state space of the problem is the 
entire /-dimensional Euclidean space 


time advance. This distance will be denoted 

U 

2. On reaching a certain distance from the run- 
way, aircraft i must assume and keep the land- 
ing speed. This distance will be denoted Z* 

These restrictions are related to the state variables 
as follows: aircraft i can no longer exercise delay or 
time advance if Si > Li , and must move at landing 
speed if Si > Z*. The latter state space construct is 
used below to formulate a class of optimal control 
problems. 

Separation constraints and landing se- 
quence 


s = r/ 

Different portions of this space will be used below 

in the optimal control models that correspond to 

different flight phases. 

The rest of this section consists of some speci- 
fies lions of I lie above machinery lo the models be- 
low. The distal point (definition 2) will be assumed 
to be a runway. Two different paths may termi- 
nate with two different, runways. The following two 
operational restrictions will be reflected: 

1. On reaching a certain distance from the run- 
way, aircraft i can no longer exercise delay or 


Denoting the minimal separation distance for an or- 
dered pair (i.i') of aircraft by p^/, and letting | • | 
denote the Euclidean norm in R 3 , we obtain the 
separation eonstraints 


|£ P (i)(siW) - z P (i')M*))l > pm' 

for all € A and all meaningful t 


( 3 ) 


If the paths of a pair (i, i') of aircraft join once and 
coincide henceforth, the set of states where these 
two aircraft are in conflict can be approximated 
from the outside by a polyhedral region, henceforth 
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Figure 2: Approximate pairwise conflict zones. 


called a pairwise, conflicting state zone and denoted 
Cij'. whose dimensions are determined by the re- 
quired minimal separation distances, and whose di- 
agonal extent is commensurate with the length of 
the common portion of the two paths. For example, 
the polygonally approximated conflict zones corre- 
sponding to the path pairs shown in Figure 1 are 
shown in Figure 2. The resulting form for the 
separation constraints is 

s ^ U ij'eA Ci# (4) 

Remark 2 The separation distance required for a 
pair of aircraft generally depends on their wake 
classes and on the order in which the two aircraft 
consecutively enter a common path. This implies 
the possibility of the asymmetry pi r 

Landing sequence and the minimal sepa- 
ration line (MSL) 

In this paper, the landing sequence is assumed pre- 
determined. (While this assumption can be justi- 
fied operationally today, it may lose basis in future 
ATM operations and research. The computational 
implications of abandoning the assumption are sum- 
marized briefly in the sect ion Discussion.) A given 
landing sequence can be described as a permutation 
7T on the set A = { 1, 2, . . . , F} such that, at the time 
when the first aircraft in the set lands and all the 


others still observe separation, the state vector s 

satisfies the inequalities 

%(1) < s tt( 2) < < s ~(/) (5) 

Thus, the permutation tt “sorts the aircraft in the 
order of reaching their distal points” , which in this 
case are the runway or runways. 

Given a landing sequence a state s with min- 
imal separation between every consecutive pair of 
aircraft satisfies the following, stronger, version of 
the inequalities (5): 

S !T(*)+/M<)-*(*+l) = ^7T(i+l)» * = 1,2,... ,7—1 (6) 

The system (6) of (7 — 1) equations in the 7- 
dimensional space R 7 determines a line, called the 
minimal separation line (MSL). 

Control law 

Ignoring inertia, the speed of aicraft i is assumed 
capable of being arbitrarily set to any value in the 
operationally admissible range. This leads to the 
following control law: 

s = V, (7) 

where 

V = («!,...,«/) 

The operationally admissible speed ranges will be 
specified for each of the control problems formu- 
lated below. 



Delay-time advance phase with 
fuel cost: optimal control problem 
PI 

Problem formulation 

Operationally, flight i is assigned a tentative sched- 
uled time of arrival (STA), denoted STA{ , which 
marks the end of the flight’s Delay-Time Advance 
phase. If this assignment were never subject to 
change, the appropriate condition for flight i to exit 
t his flight phase would be 


Si(STAi) = -Li 

However, STAs are updated in operations, typi- 
cally by requiring postponement of arrival by a time 
duration called delay, necessitated when the level 
of traffic demand reaches the throughput capacity 
and, without rescheduling, will lead to separation 
violations. Allowing not only to postpone phase 
exit, but also to expedite it (i.e., allowing a time 
advance ), carries a potential operational gain: dur- 
ing periods when traffic demand is below satura- 
tion, time advance can ultimately allow aircraft to 
be landed earlier than originally planned. Thus, 
with a corresponding adjustment and optimization 
in runway scheduling, a practice of time advance 
can help increase throughput. 

In the present model, delay and time advance 
are represented by introducing for each aircraft i a 
variable 5i whose value is to be determined as part 
of solving the optimization problem. This results in 
the more realistic condition 


SjfSTAi + di) — — Li 

for flight i to exit the Delay-Time Advance phase. 

In this paper, sufficiently early knowledge of an 
aircraft’s whereabouts (starting at a time, denoted 
by £°) is assumed. This gives the following initial 
condition: 

s (t°) = s° € S°, t° < min {STAi + Si) (8) 

i 

The quantities di are unknown at the outset and can 
be found, once solutions for all three flight phases 
are obtained, using formula (22). 


The Delay-Time Advance phase of an aircraft V s 
flight lasts until Si reaches the value —Li. Accord- 
ingly, denoting by T the exit time for the problem, 
the target set is defined for this problem as 

S r ’ p i = {s : Si < —Li for all i, with 

equality for at least one ?'.} 

For a case of 2 a/c, Figure 3 schematically shows 
the state space for this problem as the quad- 
rant { (<s*i , S 2 ) : s\ < — Li,S '2 < — L 2 } shaded 
gray, with an example of an initial state (8). 




Figure 3: A 2-a/c state space for PI. 

The speed range for aircraft i is the full cruise 
speed range , Vi\ , which contains a “economical 
cruising speed” v^ cs , closer to v than to v that 
minimizes the instantaneous rate (11) of fuel burn. 
For future reference, we write the constraints on the 
speeds: 

Vi < Vi < Vi i = 1,2,...,/ (10) 

The running cost of for aircraft Si with s* < — Li 
is the instantaneous rate of fuel burn, B r ex- 
pressed as a known function of the aircraft’s speed 
Vi. Based on [13], B r ^ is a convex function on 
\Vj,Vi] with a minimum at some “economical cruis- 
ing speed” vf cs . 

For simplicity, all vf cs ’s will be assumed to have 
the same value, henceforth denoted Q. All theoreti- 
cal results of this paper, ho wever, easily generalize to 
the case when the ’s are not necessarily equal. 

Throughout this paper, the functional form 

B r(i)i v i) = 1 + n( v i - 


( 11 ) 




with f3 = vf 28 , is used for the instantaneous rate of 
fuel burn for each i, to mimic the data shown in 
reference 1 13, figure 1|. Although symmetric about 
/3, the form (2) exhibits qualitative behavior that 
resembles the asymmetric form in [13, figure I suf- 
ficiently for our first-step investigation. 

The total running cost is, therefore, 

/ 

( 12 ) 

i= 1 

where the coefficient c measures the relative utility 
of minimizing fuel burn to minimizing excess sepa- 
ration. The functions B r (q are assumed quadratic 
polynomials, which can be justified by interpolat ing 
B T (i) through 3 data points. 

Let s MSL be the intersection between the MSL 
(6) and the target set S T ’ P1 defined in (9). To incor- 
porate the goal of minimizing separation, the exit 
cost for problem PI is defined as 


strategy, v, which is constant. (The sample nu- 
merical solutions obtained using the Dynamic Sys- 
tem Optimization Algorithm (OCP) software [14] 
are consistent with the latter conclusions.) 

Remark 3 The qualitative insight gained from the 
above application of the Pontrgagin Maximum Prin- 
ciple is that problem PI has an optimal control strat- 
egy v* which is constant. This insight not only 
shows that v* satisfies the operational preference ex- 
pressed in remark 1 , but also reduces a problem in 
optimal control to one in static optimization, 

Such a constant control strategy can be found by a 
direct polynomial- time computation as follows. 

The total cost of a trajectory s(t) corresponding 
to a control strategy 

v(t) = v = const. (15) 

is 


£(s) = i||s-s MSL || 2 , (13) 

where || • | is the norm (reference [1, section 5.2- 
5]) arising from the inner product (reference [1, 
section 5.2-6]), denoted here bv {■,■). Note that 
s t(/) L — — and the rest of the sf ISL ’s can be 
directly calculated from (6). 

Definition 2 The optimal control problem given by 
(7), (8), (9), (10), (12), (13) will henceforth be 
called problem PI. 


<T> 



T, B r(i)( v i) dt +^ | Is 1 


MSL\\2 

o 


Recalling that 

I Is 1 - s MSL || 2 = (s 1 - s AISL , s 1 - s MSL ) , 

and that, in light of ( 15 ). the final state s 1 is given 
by 

/•T 

s 1 = s° + / v dt = s° + Tv, 

Jt = o 


we obtain 


Solutions for problem PI obtained using 
the Pontryagin Maximum Principle 

Denote the costate variable corresponding to Si by 
a i and write <r = (oi, ... , 01 ). The Pontryagin func- 
tion (also known as the variational Hamiltonian ) for 
the problem is then 

H = -cJ2 B r(i)(vi) + (v. <t) (14) 


The adjoint state [3] equations 


& = -2-h = o 

as 


imply that O' is constant. From the constancy of O', 
we conclude that there exists a maximizing control 


<I> Br ( i){Vi) + jlls 0 — s M5i || 2 

+ |T 2 ||v|| 2 + T(s 0 -s MSL ,v) 

Omission of constant terms from the latter expres- 
sion leaves the solutions of the minimization prob- 
lem (min — > <h) unchanged, and therefore the 
optimal control problem can be stated as 

miii v ,r — Mh = cTYa B r(i) (vi) 


+ ^T 2 ||v|| 2 + T(s° — s MSL .v), 

(16) 


subject to the constraint 


def o 

9 = $i 


+ Tvi — (—Li) = 0 



Using a Lagrange multiplier [1] A, we arrive at the 
problem of minimizing the Lagrangian 


The function <I>. defined in (16), can be written 
<f> T a 2 (T)||vj | 2 + <a\v) + a 0 , 

where 

a 2 {T) = ac + ~T, 

a 1 = S° - S MSL - 2a ft cl, l = (l, . . . , l) e R 7 , 
and 

a f , = cF( 1 + aft 2 )- 

For each fixed value of T. the latter expression in []’s 
is a quadratic form [1] in v, whose principal axes [1] 
we now exhibit by several changes of variables, thus 
simplifying minimization. First, on defining 

u = y^rjv, 


A d 4 f cl> - A g 
by solving the algebraic system 

0 = 3§tA i = ' 

0 = 77^ A (19) 

" =9 
for w. T. A. 

For i = 2, . . . , F, one obtains 

0 = tt— A 2 Twi, 



whence Wi = 0. 

Thus, regardless of the dimension /, the above 
minimization problem reduces to the static 2- 
dimensional minimization problem 

min -> T wj - i||b 1 (7 , )|| 2 + a 0 , (20) 

wi,T 4 

subject to the constraint (18). This reduction gives 
the exact coefficients of the following polynomial in 
the variable z = V^C 77: 

- z 8 ( 16 ( 2 ( 1 ) 

+ z 7 (4L( + 4.s-? + 4aJ) 

+ z 6 (32aoac) 

+ z 5 (— 4< tel. i — 4(w.s’i) 

+ z 4 (8acl 2 - 16aoA 2 c 2 ) 

+ z 3 (8n 2 c 2 L 1 + 8a 2 c 2 s? - 12a\a 2 c ?) 

+ z(8aja 3 c 3 ) 

— 16a 2 c 2 2 2 / 2 + 8 (v 3 c 3 I 2 

The latter polynomial, obtained using the software 
[15], is readily seen to have degree 8 independently 
of 7. One of the roots of the polynomial corresponds 
to the solution of (20), (18). This gives a computa- 
tional solution procedure that has the desired fea- 
tures discussed in the latter portion of the section 
Introduction. This knowledge enables a calculation 
of the solution with arbitrary accuracy for the gen- 
eral /-aircraft problem PI in time 0(7) if the in- 
tended exit sequence for PI is pre-determined, and 
in time 0(7 2 ) otherwise. 




Obtaining delays and time advances from 
a solution 

Once a solution v(£) covering all three flight phases 
is found, the quantities 8% appearing in (8) are found 
as follows. For aircraft i, let T- 1 be the earliest time 
when 

Si = -L u 

take 

4 = T} - ST A, (22) 

Numerical experiments for prob- 
lem PI 

The correctness of the results obtained in the previ- 
ous section was tested by solving several instances 
of problem PI using the following three methods. 

1. Solving problem (20), (18) by finding numer- 
ically a suitable root of (21) and, from it, the 
corresponding speed advisory v h , using the 
wxMi axima software [15]. 

2. Solving problem (20), (18) numerically using 
the the f minsearch function in the MATLAB 
software [16]. 

3. Solving the original, /-dimensional, optimal 
control problem PI, given by (7), (8), (9), 
(10), (12), (13), using the OCP software [14] 
and without assuming the constancy of opti- 
mal control strategics asserted in Remark 3. 

The obtained values from all three approaches 
are shown in Tables 1 and 2. For each value of / 
(number of aircraft), the value of the coefficient c 
(see (12)) was varied over the values 

c = 10 n , n = -2, -1.8, -1.6 , . . . , 6 

Plots of dependence of the solution and excess sep- 
aration on the value of c are shown for 2 and 3 
aircraft in Figures 4 and 5, respectively. Figure 6 
shows MATLAB solution plots for a case of 30 air- 
craft. 


The fuel burn curve 

Throughout this section, the functional form (11) of 
instantaneous fuel burn will be used with parameter 
values 

a; = 5 x 10 -5 . = 4.5 x 1 0 2 , 

chosen to mimic the data shown in reference [13. fig- 
ure 1]. 

A case of 2 aircraft 

The numerical results in Table 1 were obtained as- 
suming that i = 1 exits the delay-time advance 
phase before the other flight does, and assuming 
the following parameter values: 

L\ = 20.0 nmi, L 2 = 30.0 nmi 

() = 5 nmi = min. separation required 
for i = 2 trailing i = 1 

> (23) 

S msl = { _ Li _ L2 _ p) 

s° = (. -L ! - 0, -L 2 - 0) 

Vj = 250 kts, 77, = 550 kts for all i 

Solution values are shown in Table 1; solution plots, 
in Figure 4. 

A case of 3 aircraft 

The numerical results in Table 2 were obtained as- 
suming that i = 1 exits the delay-time advance 
phase before the other flight s do, and assuming the 
following parameter values: 

L\ = 20.0 nmi, L2 = 30.0 nmi, 

L3 = 25.0 nmi 

f = 5 nmi = min. separation required 
for all aircraft pairs 

' (24) 

s msl = ( _ Li _ l , 2 _ A . L . s _ 2f)j 
S ° = (-Li 0,- L-, - 0 . L :i -0) 


Vj = 250 kts, Vi = 550 kts for all i 



Table 1: Solutions for instances (23) of Pi. 


c 

Method 

v* (kts) 

T (hr) 

l.e-2 

w(xMaxima) 

M(ATLAB) 

O(CP) 

(4.77e+02, 

(4.74e+02, 

(4.74e+02, 

4.71e+02) 
4.69e+02 ) 
4.69e+02) 

9.44e-01 

9.49e-01 

9.49e-01 

l.e-1 

w 

M 

O 

(4.74e+02, 

(4.74e+02, 

(4.74e+02, 

4.69e+02) 

4.69e+02) 

4.69e+02) 

9.49e-01 

9.49e-01 

9.49e-01 

l.e+0 

w 

(4.74e+02, 

4.69e+02) 

9.49e-01 


M 

(4.74e+02, 

4.69e+02) 

9.49e-01 


O 

(4.74e+02, 

4.69e+02) 

9.49e-01 

l.e+1 

w 

(4.74e+02, 

4.69e+02) 

9.49e-01 


M 

(4.74e+02, 

4.69e+02) 

9.49e-01 


O 

(4.74e+02, 

4.69e+02) 

9.49e-01 

l.e+2 

w 

(4.74e+02, 

4.69e+02) 

9.49e-01 


M 

(4.74e+02, 

4.69e+02) 

9.49e-01 


O 

(4.74e+02, 

4.69e+02) 

9.49e-01 

l.e+3 

w 

(4.75e+02, 

4.68e+02) 

9.47e-01 


M 

(4.75e+02, 

4.68e+02) 

9.47e-01 


O 

(4.75e+02, 

4.68e+02) 

9.47e-01 

l.e+4 

w 

(4.81e+02, 

4.62e+02) 

9.36e-01 


M 

(4.81e+02, 

4.62e+02) 

9.36e-01 


O 

(4.81e+02, 

4.62e+02) 

9.36e-01 

l.e+5 

w 

(4.90e+02, 

4.53e+02) 

9.19e-01 


M 

(4.90e+02, 

4.53e+02) 

9.19e-01 


O 

(4.90e+02, 

4.53e+02) 

9.19e-01 


Solution values are shown in Table 2; solution plots, 
in Figure 5. 

A case of 30 aircraft 

The timeframe in which the 2- and 3-aircraft in- 
stances of PI were solved, above, is on the order 
of 1 hour. Within this timeframe, saturation rates 
of air traffic arriving in a terminal can well exceed 
60 aircraft per hour, while minimal rates are typ- 
ically far in excess of 3 aircraft. Thus, a case of 
30 aircraft corresponds to a rate of arrivals which 
is, on the one hand, below saturation and, on the 
other hand, more realistic than 3 aircraft. It is at 
this rate that the advantages of time advance are 
expected to be most manifest. 

The numerical results for this case were obtained 
assuming that i = 1 exits the delay-time advance 
phase before the other flights do, and assuming the 


following parameter values: 

Lj = 20.0(1 + t<i) nmi, where m is 
a random variable [1] distributed 
uniformly [1] in [—0.3, 0.3] 

p = 5 nmi = min. separation required 
for all aircraft pairs 

► 

s MSL = (sf* SL ) h sf ISL = - Li - (i - 1 )ft 

s° = (—Li fi)i for all i 

v_i = 250 kts, Vi = 570 kts for all i 

The speed advisories were computed using 
MATLAB [16] without imposing the speed range 
constraints. For values of c above 10 4 , the com- 
puted speed advisory u* for aircraft i = 1 exceeded 
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Figure 4: Solutions for the 2 a/c case of PI. 


the value v * = 570 kts. This behavior may be an 
indication that the functional form B r ^(ui) of fuel 
burn used above does not include the highly nonlin- 
ear effects arising as a subsonic aircraft approaches 
the speed of sound, and also that ascribing exces- 
sively large priority c to fuel savings can yield op- 
erationally infeasible results. This issue will be ad- 
dressed as the model and algorithm are further de- 
veloped to accommodate a higher level of opera- 
tional realism. For c < 10 4 , the computed speed 
advisories do meet the speed range constraints, the 
values of E(s) are below 5 nmi, and fuel burn ex- 
hibits a negative correlation with excess separation 
(see Figure 6). Consequently, fuel savings are neg- 
atively correlated with throughput, as expected. 


Discussion 

Summary: the results and their opera- 
tional implications 

The optimal control formulated above (problem PI) 
describes aircraft t raffic in its Delay - Time Advance 
phase of traffic arrival. The cost functional includes 
two terms, collective for the entire traffic: fuel burn 
(12) and excess separation (13). These two terms 
are in competition for priority; e.g., in order to re- 
duce excess separation (thereby increasing through- 
put), some aircraft must speed up, incurring higher 
fuel burn. The cost functional thus represents the 
tradeoff between the pursuits of minimizing the two 
costs. The relative weight ascribed to fuel burn over 
excess separation is represented by the coefficient c 
appearing in (12). 

Similar optimal control models are formulated 




Table 2: Solutions for instances (24) of Pi. 


c 

Method 

v* (kts) 

T (hr) 

l.e-2 

w(xMaxima) 

M(ATLAB) 

O(CP) 

(1.77e+02. 

(4.77e+02, 

(4.77e+02, 

4.72e+02, 

4.72e+02, 

4.72e+02, 

4.66e+02) 

4.66e+02) 

4.66e+02) 

9.43e-01 

9.43e-01 

9.43e-01 

l.e-1 

w 

M 

0 

(4.77e+02, 

(4.77e+02, 

(4.77e+02, 

4.72e+02, 

4.72e+02, 

4.72e+02, 

4.66e+02) 

4.66e+02) 

4.66e+02) 

9.43e-01 

9.43e-01 

9.43e-01 

l.e+O 

w 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 


M 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 


0 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 

l.e+1 

w 

(4.77e+02, 

4.72e+02, 

4.66e+02, 

9.43e-01 


M 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 


0 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 

l.e+2 

w 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 


M 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 


0 

(4.77e+02, 

4.72e+02, 

4.66e+02) 

9.43e-01 

l.e+3 

w 

(4.78e+02, 

4.71e+02, 

4.66e+02) 

9.41e-01 


M 

(4.78e+02, 

4.71e+02, 

4.66e+02) 

9.41e-01 


0 

(4.78e+02, 

4.71e+02, 

4.66e+02) 

9.41e-01 

l.e+4 

w 

(4.87e+02, 

4.65e+02, 

4.63e+02) 

9.24e-01 


M 

(4.87e+02, 

4.65e+02, 

4.63e+02) 

9.24e-01 


0 

(4.87e+02, 

4.65e+02, 

4.63e+02) 

9.24e-01 

l.e+5 

w 

(5.06e+02, 

4.54e+02, 

4.54e+02) 

8.90e-01 


M 

(5.06e+02, 

4.54e+02, 

4.54e+02) 

8.90e-01 


0 

(5.06e+02, 

4.54e+02, 

4.54e+02) 

8.90e-01 


(see appendices I and II) to describe the two sub- 
sequent phases of each arrival flight. The solution 
in each model is a collective speed v(£), which is 
generally a time-dependent vector whose i - th com- 
ponent is the speed advised for flight i along its 
intended path. (As is shown above, v(£) turns out 
constant for the first, Delay - Time Advance, phase 
of arrival.) Once the solutions v(t) for all three 
flight phases are found, t hey completely determine 
the optimal (in the sense of the cost functionals) 
amount of delay or time advance for each aircraft. 
This determination is given by formula (22). 

The solutions obtained for the Delay - Time Ad- 
vance phase (problem PI) and verified using three 
independent computational methods (tables 1 and 
2) are a step toward automation of scheduling that 
takes account of fuel burn and its interplay with 
excess separation. The latter quantity, in turn, af- 


fects airport efficiency. The plots shown in Figures 
4 and 5 indicate that, i)as expected, with the higher 
“relative importance” (value of c) ascribed to fuel 
burn, excess separation increases toward, and levels 
out at, a maximal value, and ii) the speed advised 
for the lead aircraft increases. The plot in Figure 
6 confirms the intuitive expectation that increasing 
the priority of fuel savings results in higher excess 
separation, hence in lower throughput. 

Another potential benefit, of such quickly and 
accurately computable solutions is that, with 
equipage to communicate the parameters q, ft. Li, s° 
to the crew of each flight i and to carry out simple 
calculations, each flight crew is in a position to de- 
termine its optimal speed advisory that meets the 
separation constraints between this flight and all 
the others. The contributions of this paper, conse- 
quently, are a step toward self- separation [12]. 




</) 


Optimal speeds vs. log 1Q c. 
Legend: ’o’ v^, V v*, ’+’ 


500 h 


’> 480^- 

’ >T 460 
II 

' > 


ooooo 


6000000000 ( 1 ) 


O O O D O D O o o o O O-O-O -O. ... 
******* $$$ ^ $ ^^ 


++++++ 


* * 




0 


lo 9io c 

log 1Q (excess separation) vs. log in c 


'10 


Uj 


S 5 


2 r 


Oh 


-2 





.. A . .< 

0 0 

0 0 

> 0 0 

V 

o 

o 


...... a. .. 5 

o 

o 

. . .A 

0 0 





o 

o 

0 v 










2.5 


3.5 


lo 9io c 


4.5 


5.5 


^ 115 
^ 1.1 

e 105 

03 c ! 

U) O) 0.95 

> ^ 

03 ^ 


Avg fuel burn vs. (excess separation ) 1 


. 1/2 











=1 :e06 

CO... _ 

c-: i. exOO. . . 






O O- 

c 

■ Q.. O.c 

<KOO O 

o c 

o 

o 

6 

c 










1 1 1 1 1 1 1 1 1 1 


10 15 20 25 30 35 40 45 

E( s ) 1/2 (nmi) 


Figure 5: Solutions for the 3 a/c case of PI. 


Directions for future research 

The following research directions seem sensible for 
next steps toward operational realism and develop- 
ment of automation tools lit for use in the field. 
Introducing uncertainty terms into the control laws 
describing the traffic can help capture the behavior 
of traffic in the presence of wind and human fac- 
tors in control execution. A solution the 3-phase 
control problem with PI, P2, and P3 as phases 
would enable formula (22) to be used for calculat- 
ing the resulting delays and time advances. Inclu- 
sion of inertia in the control model can indicate how 
the cost function may need to be modified in order 
to capture the resources of importance (e.g., fuel) 
and to yield speed advisories that, like the ones 
found above, are piecewise constant (Remark 1). 
Testing the computational procedures using high- 
fidelity ATM operations and live traffic. 
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Appendix I: The optimal control 
problem for the Mixed Phase 

This problem, denoted P2, arises once the system 
reaches a state s 1 in the target set (9) of problem 
PI. Given the state s 1 , which will serve as the ini- 
tial state for problem P2, let the set A\ consist of 
those indices i in A for which the aircraft still has 
the option of exercising delay or time advance, i.e. 

«! < -Li 

Problem P2, while governed by the same control 
law as PI, allows aircraft i a narrower speed range 

The target set for P2 is the set of all states s for 
which 

Si = —Li for at least one index i in Ai 

The running cost includes a term reflecting fuel 
burn only for those aircraft in A\, but the exit cost 
is the same as in PI. The resulting Pontryagin func- 
tion is 

H = -c D r(i)( v i) + <v-cr) 
feAi 

Appendix II: The optimal control 
problem for the Separation Mini- 
mization Phase 

Once each Si exceeds —Li, none of the aircraft pur- 
sues the objective of minimizing fuel burn and is 
concerned solely with minimizing separation. The 
control, the speed ranges, and the exit cost are as 
in P2. The running cost is now zero. The target set 
is a state s 3 in which 

—Li < Si < —li 

The Pontryagin function is 

H = (v, cr) 
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